Energetics of the oxidation and opening of a carbon nanotube 
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We apply first principles calculations to study the opening of single-wall carbon nanotubes 
(SWNT's) by oxidation. We show that an oxygen rim can stabilize the edge of the open tube. 
The sublimation of CO2 molecules from the rim with the subsequent closing of the tube changes 
from endothermic to exothermic as the tube radius increases, within the range of experimental fea- 
sible radii. We also obtain the energies for opening the tube at the cap and at the wall, the latter 
being significantly less favorable. 



PACS Numbers: 61.48. +c, 71.20.Tx, 71.15.Nc 

The possibility of filling carbon nanotubes jjj with a 
variety of substances has fascinated the materials science 
community since they were first observed Q and synthe- 
sized in promising quantities [0], but the fact that the 
tube tips are closed represented a severe limitation for 
it. Ajayan and Iijima Q observed the cap removal and 
filling of the nanotubes when they were heated in air, 
suggesting a chemical reaction between oxygen and the 
tube ends. This was confirmed by Ajayan et al || who 
showed that the tubes were open exclusively at the caps. 
Simultaneously, Tsang, Harris and Green described 
the thinning and opening of nanotubes by carbon diox- 
ide. They also observed corrosion at the cap region, the 
tube walls remaining intact. Ugarte et al JtJ described a 
method of filling nanotubes in which they were first ox- 
idized in air to remove the caps. In this process, 60% of 
the tubes were open. The above mentioned experiments 
used multi-walled nanotubes. In a recent experiment, 
however, Dillon et al. [|j observed the physisorption of H2 
in bundles of small-diameter SWNT's. They suggested 
that these tubes were open by oxidation in the same way 
as multi-walled tubes. Individual open-ended SWNT's 
were also obtained experimentally || by treating SWNT 
rope material with a mixture of acids. 

The filling and opening of nanotubes have also been 
subjects of theoretical studies. Pederson and Broughton 
fiof predicted the tendency of the nanotubes to attract 
molecules to their interior. More recently, Miyamoto et 
al jll[] reported first-principles calculations that showed 
the ionic cohesion between carbon nanotubes and potas- 
sium atoms linearly arranged inside them. The reactivity 
of the open end of SWNT's was studied by Charlier et 
al jL2|. They showed that the unsaturated tip of small 
diameter SWNT's closes spontaneously at typical growth 
temperatures. However, first principles calculations also 
indicate that these tubes may stay open during growth 
in the presence of a catalyst, such as Ni JT3| . 

The detailed description of the burning and opening 
processes, and of their kinetics, is a formidable problem 
beyond the scope of state-of-art ab initio techniques. In 



the present work, we apply first principles calculations to 
address two specific, but relevant questions on the ener- 
getics of these processes: (i) why does O2 attack pref- 
erentially the tube caps? and (ii) how and why do the 
tubes remain open (if at all) after the O2 supply has been 
cut? To address the first question, we compare the open- 
ing of oxidized holes at the cap of the tube with their 
chemical equivalents at the wall, finding that the cap ox- 
idation is much more favorable, due to the release of its 
curvature strain. On the second question, we find that 
the presence of an oxygen rim at the edge of the open 
tube results in a very stable configuration, even in the 
absence of 02- Smaller holes with favorable energetics 
allow us to obtain an estimate of the activation barrier 
of the closing process. 

Our calculations were performed with the SIESTA 
program ]I^16|, using density- functional theory Jl7| , 
within the generalized gradient approximation (GGA) 
for exchange-correlation p8| , and norm-conserving pseu- 
dopotentials Jn|^o|]. The basis set is a linear combina- 
tion of pscudoatomic orbitals |2l| , |2^|. We always used a 
split-valence double-^ (DZ) basis |2j|, and convergence 
tests were performed by adding polarization functions 
p4| . This resulted in changes in formation energies of 
only a few meV/atom. All the geometries were opti- 
mized, with remanent forces of less than 0.1 eV/A. Test 
calculations were performed on graphite, nanotubes, and 
fullerenes. The calculated nearest-neighbor distances of 
C 6 o (1-42 and 1.47 A) and graphite (1.42 A) reproduce 
the experimental values (1.40, 1.46, and 1.42 A, respec- 
tively) ||] within 2%. Also, the energy per atom of Ceo 
relative to that of graphite (0.38 eV/atom) agrees very 
well with the experimental result (0.392 - 0.433 eV/atom) 
p5| . Our calculations were performed for (4,4), (6,6), and 
(8,8) armchair SWNT's, with diameters 5.6, 8.3, and 11.1 
A, containing up to 166 carbon atoms. However, for the 
sake of clarity in the exposition, we will refer to the spe- 
cific (6,6) case throughout the paper, leaving the analysis 
of the size effects to the end. 

In order to understand why O2 produces a selective ox- 
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FIG. 1. Structures used in the calculations for the (6,6) 
carbon nanotube. (a) Infinite tube (periodic boundary con- 
ditions); (b) closed tube, with 120 carbon atoms and 12 hy- 
drogen atoms saturating the bulk-like end of the tube; (c) 
oxidized open tube; (d) oxidized opening equivalent to (c) 
but at the wall instead of the cap; (e) oxidized hole obtained 
by removing an hexagon at the cap; (f) same as (e) at the 
wall; (g) oxidized hole obtained by removing a pentagon; (h) 
same as (g) but removing an additional C atom; (i) a C va- 
cancy at the cap; (j) vacancy at the wall; (k) unsaturated 
open tube; (1) H-saturated open tube on both sides (used to 
characterize the saturation energy). The free Hb, O2, H2O, 
and CO2 molecules were also calculated. The back atoms are 
not shown in the figure. 

idation of the tube caps, rather than a generalized burn- 
ing, we compare the energies of similar holes at the tube's 
cap and wall. Figure 1 shows these structures for the 
(6,6) tube. The closed tube, hereafter denoted by C120, 
is shown in Fig. 1(b). The other end of the tubes is bulk- 
like and saturated by hydrogen atoms. This termination 
is used in all the structures, except the infinite tubes used 
to calculate the carbon chemical potential \xc-, for which 
we used periodic boundary conditions [Fig. 1(a)]. The 
hole geometries, shown in Fig. l(c)-(j), were chosen as 
plausible low energy strutures for specific hole sizes and 
not necessarily to simulate the actual oxidation process. 



These hole structures will provide estimates of the en- 
ergy released, for a given hole size, during the tube open- 
ing process. Figure 2 shows the potential energy change 
during oxidation as a function of the hole size, measured 
by the number of carbon atoms lost by the tube in the 
process. These atoms are removed from the tube accord- 
ing to the reaction C+O2 — > CO2, while extra oxygen is 
used to saturate the resulting holes. The initial structure 
is Fig. 1(b) and the oxidized structures are Fig. l(c)-(j). 
The general trend is linear with the number of C atoms 
missing in the hole, as it should, but the energy released 
by oxidizing the cap is larger than for the wall, as shown 
in Fig. 2(b) ^6|. As we will see, the origin of the differ- 
ence is the extra elastic energy accumulated at the cap. 
It is remarkable that the wall-cap difference for remov- 
ing the first carbon atom (2.4 eV), is already very close 
to that for a full hexagonal hole. Thus, the cap's strain 
energy is substantially released from the very beginning 
of the oxidation. 

Addressing the stability of the open tube after oxida- 
tion, Fig 1(c) shows the most efficient chemical saturation 
of the open tube with oxygen. In this oxidized open tube, 
denoted by C120O6, every O atom saturates two carbon 
"dangling bonds", giving a slightly stretched C-0 bond 
distance of 1.43 A and a C-O-C angle of 105°. The bond 
stretching produces an inward relaxation (18°) of the rim 
structure as a consequence of the pulling of the O atoms, 
diminishing the tube perimeter at the rim. The struc- 
ture is very stable: the sublimation of oxygen atoms in 
the form of CO2, resulting in an unsaturated open tube 
[Fig. l(k)], is found to be energetically unfavorable by 
17.8 eV. This highly endothermic reaction indicates that, 
if open tubes are present, they will have an oxygen rim 
at the edge even at typical oxidation temperatures (~ 
1000K) fg. 

Starting with this oxidized open tube, we consider a re- 
action in which the tube closes and CO2 molecules evap- 
orate: 

oxidized open tube — > closed tube + 3CO2. (1) 

The energy of this reaction is AE = (Ec 120 + 3i?co 2 ) — 
( E C iao o e + 3^c) The result is AE = 0.2 eV with the DZ 
basis, and AE — 0.3 eV including polarization orbitals. 
Thus, the oxidized open tube is slightly favorable ener- 
getically over the closed one for this tube diameter. How- 
ever, it has to be stressed that the system composed by 
the closed tube and the carbon dioxide gas has larger en- 
tropy and will thus predominate in thermodynamic equi- 
librium at a high enough temperature. 

In addition to its genuine stability, the oxidized open 
tube can be further stabilized by the energy barrier of 
the closing process, which requires substantial rebond- 
ing. The actual rebonding is extremely hard to guess 
or simulate. We search for an estimate of the energy 
barrier as follows: first we consider a "closing coordi- 
nate" that the system must necessarily follow during the 
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FIG. 2. (a) Potential energy change, E;, urn , during the 
burning process of the original closed tube, leaving behind 
the holes shown in Fig. 1. Ef, urn is plotted versus the number 
of C atoms missing in the final hole. Filled symbols are for 
holes at the cap of the tube, open symbols for holes at the 
wall. The labels refer to the structures of Fig. 1. (b) Energy 
difference between burning the tube at the wall or at the cap: 

Ei rpwall r^cap 



closing process. We simply use the hole size, denned as 
before as the number of carbon atoms missing from the 
closed tube. Then, for particular values of this coor- 
dinate, we search for low energy structures using basic 
chemical principles: to break the minimum number of 
the most strained C-C bonds, and to saturate them with 
oxygen. The chosen structures are those in Fig. 1(e), (g), 
and (h), obtained by removing an hexagon, a pentagon, 
and a pentagon plus an extra C atom, respectively, and 
introducing three O atoms saturating the holes. We shall 
mention that other small-hole structures or other closure 
mechanisms (with catalysts, for instance) with even lower 
energies can possibly exist, our construction giving just 
a plausible estimate of a low energy barrier. Again, we 
considered reactions in which the small holes close and 
CO2 molecules evaporate, and we evaluated the energy 
difference between the initial and final states. The re- 
sults are summarized in Fig. 3, which shows (solid line) 
these hole energies as a function of the number of car- 
bon atoms missing in the hole. We can see in the figure 
that for closing the tube just 0.2 eV are needed, but a 
barrier of at least 4.2 eV has to be overcome first. Re- 
markably, the hole formation energies for the three small 
holes considered here are degenerate within 0.1 eV. As 
for the equivalent wall holes, displayed in Fig. 1(d) and 
(f), Fig. 3 shows (dashed line) that the tube closing is 
clearly exothermic (by 9.2 eV) and that no energy bar- 
rier is found. 



We can understand the origin of the distinct cap and 
wall behaviors by defining a cap energy E cap , as the en- 
ergy of the capped tube minus that of the same number 
of atoms in the bulk of the tube. Since chemical bonds 
are similarly saturated at the cap and at the wall, the ori- 
gin of E cap must be the larger elastic energy of the cap. 
It amounts to 9.0 eV for the (6,6) tube, which is com- 
parable to the 9.2 eV energy obtained for the opening of 
the large hole at the tube wall. 

Since the release of the strain energy, which plays a 
crucial role in the opening process, is related to the tube 
diameter, narrower tubes should be more easily open due 
to their larger cap curvature. Our calculations on the 
(4,4) and (8,8) armchair tubes confirm this trend: with 
CO2 molecules used as an oxygen atom reservoir, as in 
equation (1), the closing energy amounts to 1.2 eV for 
the (4,4) tube, 0.2 eV for the (6,6) tube (as we have 
already shown), and -0.3 eV for the (8,8) tube, allowing 
us to estimate a value of ~ —0.6 eV for the (10,10) tube. 
Therefore, within the range of experimental feasible radii 
p8[ the CO2 sublimation changes from endothermic to 
exothermic, always within tenths of eV. 

The presence of an oxygen rim may have a strong in- 
fluence on the ability to fill the tubes. In some cases, it 
may be desirable to replace it by a possibly more inert 
hydrogen rim. Therefore, we have studied the energetics 
of reducing the oxidized open tube with H2. Simply re- 
placing O with H, followed by O2 evaporation, is slightly 
endothermic (by 0.3 eV for the (6,6) tube). However, the 
reaction becomes strongly exothermic, by 9.1 and 9.5 eV, 
when we consider CO2 and H 2 as byproducts, respec- 
tively. Thus, we expect that the reduction of the oxygen 
rim should be quite possible. We shall mention that the 
saturation of the tube ends with other molecular termi- 
nations rather than O or H could result in even lower 
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FIG. 3. Energy differences between initial and final states 
of the closing reaction: tube with hole — > closed tube + CO2. 
The energies are shown as a function of the number of carbon 
atoms missing in the hole. Solid and dashed lines are used 
for holes at the cap and wall, respectively. The labels refer to 
the structures of Fig. 1. 
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formation energies. 

In conclusion, we have studied the energetics of open- 
ing SWNT's by oxidation. We show that an oxygen rim 
can stabilize the edge of the open tube. The sublimation 
of CO2 molecules from the rim with the subsequent clos- 
ing of the tube changes from endothermic to exothermic 
as the tube radius increases, within the range of exper- 
imental feasible radii. Our results also show that the 
opening reaction occurs preferably at the tube cap, due 
to the release of its strain energy. 



ACKNOWLEDGMENTS 



We acknowledge support from the Brazilian agencies 
CNPq, FAPEMIG, CAPES, and PRONEX-MCT, and 
from Spanish DGES grant PB95-0202. P.O. is the recip- 
ient of a Motorola PCRL Sponsored Research Project. 
We also thank R. W. Nunes for a critical reading of the 
manuscript. 



[1] M. Dresselhaus, G. Dresselhaus, and P. C. Eklund. Sci- 
ence of Fullerenes and Carbon Nanotubes (Academic 
Press, San Diego, 1996). 

[2] S. Iijima. Nature. 354, 56 (1991). 

[3] T. W. Ebbesen, P. M. Ajayaxi.Nature. 358, 220 (1992). 

[4] P. M. Ajayan, S. Iijima. Nature. 361, 333 (1993). 

[5] P. M. Ajayan, T. W. Ebbesen, T. Ichihashi, S. Iijima, K. 

Tanigaki, H. Hiura. Nature. 362, 522 (1993). 
[6] S. C. Tsang, P. J. F. Harris, M. L. H. Green. Nature. 

362, 520 (1993). 
[7] D. Ugarte, A. Chatelain, W. A. de Heer. Science. 274, 

1897 (1996). 

[8] A. C. Dillon, K. M. Jones, T. A. Bekkedahl, C. H. Kiang, 

D. S. Bethune, M. J. Heben. Nature. 386, 377 (1997) 
[9] J.Liu, A.G Rinzler, H. Dai, J.H. Hafner, R.K. Bradley, 
P.J. Boul, A. Lu, T. Iverson, K. Shelimov, C.B. Huffman, 
F. Rodriguez-Macias, Y. Shon, T.R. Lee, D.T. Colbert, 
R. E. Smalley. Science. 280, 1253 (1998). 

[10] M. R. Pederson and J. Q. Broughton. Phys. Rev. Lett. 
69, 2689 (1992). 

[11] Y. Miyamoto, A. Rubio, X. Blase, M. L. Cohen, S. G. 
Louie. Phys. Rev. Lett. 74, 2993 (1995). 

[12] J. Charlier, A. De Vita, X. Blase, R. Car. Science 275, 
646 (1997). 

[13] Y. H. Lee, S. G. Kim, D. Tomanek. Phys. Rev. Lett. 78, 
2393 (1997). 

[14] The SIESTA program is described in Ref. (h|. It builds 
the DFT selfconsistent Hamiltonian in order-iV opera- 
tions. The Hamiltonian problem can be solved either with 
diagonalization (order- iV 3 ) or with order- N algorithms 
according to convenience. In this study both options were 
used. 



[15] P. Ordejon, E. Artacho, J. M. Soler. Phys. Rev. B. 53, 
R10441 (1996). 

[16] D. Sanchez-Portal, P. Ordejon, E. Artacho, J. M. Soler. 

Int. J. of Quantum Chem.. 65, 453 (1997). 
[17] W. Kohn, L. J. Sham. Phys. Rev. 140, A1133 (1965) 
[18] J. P. Perdew, K. Burke, M. Ernzerhof. Phys. Rev. Lett. 

77, 3865 (1996). 
[19] N. Troullier, J. L. Martins. Phys. Rev. B. 43, 1993 (1991). 
[20] L. Kleinman, D. M. Bylander. Phys. Rev. Lett.. 48, 1425 

(1982). 

[21] O. F. Sankey, D. J. Niklewski Phys. Rev. B. 40, 3979 
(1989). 

[22] The range of the orbitals is defined in a balanced way by a 
common orbital-confinement energy of 0.14 eV VM. The 
resulting radii are 4.52 and 5.52 for the s and p orbitals 
of C, respectively, 5.20 for the s orbital of H, and 3.56 
and 4.35 for s and p orbitals of O, all of them in atomic 
units. 

[23] The split-off function reproduces the tail of the original 
orbital beyond some radius, matching to a smooth r l (a — 
br 2 ) towards the nucleus. The matching radius is chosen 
so that the split orbital carries 15 % of the norm of the 
original one [24|. 

[24] D. Sanchez-Portal, P. Ordejon, E. Artacho, J. M. Soler, 
to be published. 

[25] H.S. Chen, A.R. Kortan, R.C. Haddon, M. L. Kaplan, C. 
H. Chen, M. Mujsce, H. Chou, and D. A. Fleming. Appl. 
Phys. Lett. 59, 2956 (1991). 

[26] Note that the numbers in Fig. 2(a) depend on the pro- 
posed oxidation process, and that other possibilities, like 
involving C and O interstitials are plausible, and would 
affect the results. The origin of the ambiguity is the pres- 
ence of three references (bulk C, CO2, and O2) for defin- 
ing the chemical potential of just two species, C and O. 
The process proposed involves no extra C atoms. How- 
ever, the energy differences shown in Fig. 2(b) are inde- 
pendent of these assumptions. There is no ambiguity in 
the case of the closing process, since, in the absence of 
O2, both chemical potentials are well defined. 

[27] At growth temperatures (~ 3000K) , molecular dynam- 
ics simulations have shown that narrow SWNT's close 
spontaneously [ fl2| . 

[28] S. Bandow, S. Asaka, Y. Saito, A. M. Rao, L. Grego- 
rian, E. Richter, P. C. Eklund, Phys. Rev. Lett. 80, 3779 
(1998). 



4 



